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Abstract 

A historically unprecedented mountain pine beetle (MPB) outbreak affected 
western Montana during the past decade. We examined radial growth rates (AD 
1860-2007/8) of co-occurring mature healthy and MPB-infected ponderosa pine 
trees collected at two sites (Cabin Gulch and Kitchen Gulch) in western Montana 
and: (1) compared basal area increment (BAI) values within populations and 
between sites; (2) used carbon isotope analysis to calculate intrinsic water-use 
efficiency (iWUE) at Cabin Gulch; and (3) compared climate-growth responses 
using a suite of monthly climatic variables. BAI values within populations and 
between sites were similar until the last 20-30 years, at which point the visually 
healthy populations had consistently higher BAI values (22-34%) than the MPB- 
infected trees. These results suggest that growth rates two-three decades prior to 
the current outbreak diverged between our selected populations, with the 
slower- growing trees being more vulnerable to beetle infestation. Both samples 
from Cabin Gulch experienced upward trends in iWUE, with significant regime 
shifts toward higher iWUE beginning in 1955-59 for the visually healthy trees 
and 1960-64 for the MPB-infected trees. Drought tolerance also varied between 
the two populations with the visually healthy trees having higher growth rates 
than MPB-infected trees prior to infection during a multi-decadal period of 
drying summertime conditions. Intrinsic water-use efficiency significantly 
increased for both populations during the past 150 years, but there were no sig- 
nificant differences between the visually healthy and MPB-infected chronologies. 



Introduction 

The ongoing mountain pine beetle (MPB; Dendroctonus 
ponderosae Hopkins) epidemic affecting pine species in 
western North America is the largest in recorded history 
(Mitton and Ferrenberg 2012). During 2001-2010, over 
2.3 million hectares of forests in Montana were affected 
by the MPB (Chaney 2011), which caused mortality to 
numerous species of the Pinus genus including lodgepole 
{P. contorta), ponderosa (P. ponderosa), and whitebark 
{P. albacaulis) pines (Logan et al. 2010). Mortality rates 
peaked in 2009 when approximately 1.5 million hectares 
of pines experienced MPB-caused death. The MPB is 
native to western North America, with occasional years 
marked by outbreak populations (Brunelle et al. 2008) 
triggered by favorable periodic climatic and/or stand 
conditions (Gibson et al. 2009), but typically only low 



levels of infestation and subsequent tree mortality occur 
in a given year. 

The magnitude of the current epidemic has been attrib- 
uted to multi-decadal climate change (Logan and Powell 
2001; Akuma et al. 2008; Bentz et al. 2010) directly 
through a reduction in low temperature-induced mortal- 
ity caused by extended cold-air incursions, particularly 
during the autumn and spring, and indirectly through 
increased drought frequency. Carbon assimilation is 
impacted by reduced stomatal conductance and reduced 
photosynthate allocations for defense mechanisms or tis- 
sue repair during drought periods, making the infected 
trees more susceptible to mortality. Additionally, warmer 
conditions have allowed a change from univoltine to biv- 
oltine life cycles (one to two generations) in some broods 
of MPB (Mitton and Ferrenberg 2012). Extensive fire- 
suppression activities that decrease the number and extent 
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of stand-replacing fires can also affect forest susceptibility 
to MPB outbreaks (Kulakowski et al. 2012). 

During the summer months, MPB attack trees en-mass 
and girdle the phloem that in turn decreases photosyn- 
thate transport within a tree and reduces tree health (Gib- 
son et al. 2009). MPB also carry the spores of two "blue 
stain" fungal species, Ophiostoma montium and Grosman- 
nia clavigera, that are inoculated into the trees by adult 
beetles. The fungi subsequently spread in the sapwood 
and phloem of infected trees, reducing water transport 
and altering resin flow in a fashion that contributes to an 
increase in tree mortality (Gibson et al. 2009). Mass- 
attacked trees typically die the following growing season, 
although trees subject to partial or "strip" attacks may 
survive several years post infestation (Gibson et al. 2009). 

Despite the widespread mortality in forest stands, not 
all trees are affected, and the reasons for variations in 
host susceptibility remain uncertain. Tree vigor, a mea- 
sure of wood production/leaf area, is significantly inver- 
sely correlated with individual tree susceptibility (Larsson 
et al. 1983; Mitchell et al. 1983) and is influenced by 
stand density, pathogenic agents, and fluctuations in pre- 
cipitation and temperature (Gibson et al. 2009). Ponder- 
osa pine trees are most likely to be susceptible in dense 
stands (>5.6 m 2 basal area/hectare) of adult tress 
(>80 years) with diameters of 20-30 cm (Gibson et al. 
2009). Tree mortality may increase consistently up to 
diameters of approximately 23 cm, with mortality of lar- 
ger trees being more random (McCambridge et al. 1982). 

Here, we compare radial growth rates of co-occurring 
populations of mature healthy and MPB-infected ponder- 
osa pine trees collected at two sites in western Montana. 
Our objectives were to: (1) compare basal area increment 
values within populations and between sites; (2) examine 
the potential influence of changes in intrinsic water-use 
efficiency; (3) compare climate-growth responses; and (4) 
postulate on what increases individual tree susceptibility 
to MPB-induced mortality within a stand where other 
growth-influencing factors are held constant. 

Methods 

Tree-ring data collection 

We collected samples in western Montana at Cabin Gulch 
Research Natural Area 28 km NE of Helena and Kitchen 



Gulch 32 km SE of Missoula (Table 1). The two sites 
were selected because they had a history of minimal 
anthropogenic impacts including livestock grazing, wood 
cutting, and fire suppression (pers. comm., Steve Shelly, 
USFS Rl RNA Coordinator) and contained both visually 
healthy and MPB-infected trees of multiple age classes 
growing in park-like stands (Fig. 1). At each site, we sam- 
pled visually healthy and MPB-infected trees using stan- 
dard dendroecological field techniques, selecting older 
trees (>150 years) at each site using visual old-age charac- 
teristics. Only open-grown trees (Fig. 1) were cored to 
ensure that radial growth was not affected by canopy clo- 
sure during, at least, the past century. Both sites are 
located on steep (25-40%), south-facing slopes with thin 
soils ranging from 20 to 50 cm (authors' observations). 

We identified visually healthy trees by the presence of 
green needles and no evidence of pitch tubes, boring dust, 
and live MPB eggs, pupae, larvae, or adults. MPB-infected 
trees were selected by the presence of red foliage in the 
crown and upon coring, blue-stained sapwood caused by 
the blue stain fungi. We additionally examined beetle gal- 
lery patterns on selected trees to confirm MPB activity 
opposed to other boring beetles. For each tree, we col- 
lected two core samples using increment borers at 
approximately 1.4-m height. Additionally, we measured 
crown height and basal diameter for each tree. 

We processed collected cores using standard dendro- 
chronology procedures (Stokes and Smiley 1968). We 
cross-dated the cores using the list method (Yamaguchi 
1991), and confirmed dating accuracy using the diagnos- 
tic program COFECHA (Holmes 1983). We used raw 
ring-width data to calculate basal area increment (BAI, 
cm 2 yr _1 ) for all cores from 1860 to 2007/8 for each of 
the four chronologies (CGR [visually healthy and devel- 
oped from 12 trees and 20 cores] and CGI [MPB- 
infected; 10 trees, 18 cores] for Cabin Gulch and KGR 
(21 trees, 39 cores), and KGI (16 trees, 30 cores) for 
Kitchen Gulch). To account for the potential influences 
of variations in tree size and age between chronologies 
from the same site, we correlated unstandardized BAI val- 
ues with standardized BAI values for each chronology. 
We found negligible differences between the chronologies 
(all r-values > 0.97) and then retained unstandardized 
BAI values for our remaining analyses as this allowed for 
direct comparisons of BAI values within and between our 
study sites. 



Table 1. Site characteristics for Cabin Gulch and Kitchen Gulch. 



Site Name 


Latitude (°N) 


Longitude (°VV) 


Elevation (m) 


Climate division 


DBH- visually healthy (cm) 


DBH-infected (cm) 


Cabin Gulch 


47.79 


111.77 


1140 


Montana 4 


32.44 


27. S3 


Kitchen Gulch 


46.71 


113.65 


1444 


Montana 1 


34.4 


35.95 
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Figure 1. Visually healthy and MPB-infected trees at Kitchen Gulch, 
Montana. Photograph by authors, August 201 1. 

Carbon isotope analysis 

We developed carbon isotope chronologies from both the 
visually healthy and MPB-infected trees at Cabin Gulch 
following methods discussed in Knapp and Soule' (2011). 
For each chronology, we selected six trees based on the 
longest and clearest ring structure, with two samples rep- 
resenting each tree. We manually separated tree rings for 
each sample into pentads (5-year segments) beginning in 
AD 1800 or later (because of lack of older cores) using a 
scalpel. We pooled the pentads from all samples and 
ground together to 40 mesh. To quantitate inter-tree iso- 
topic variability, we also pooled samples for each 50-year 
period. We processed all ground samples to a-cellulose 
following conversion to holocellulose using the Jayme- 
Wise method (see Leavitt and Danzer 1993) and using 
batch processing and commercial digestion pouches (AN- 
KOM Technology, Boston, MA). Sample extractions were 
first removed using a soxhlet extraction apparatus operat- 



ing with toluene/ethanol and ethanol organic solvents fol- 
lowed by boiling in deionized (DI) water. Lignin was 
removed by reaction in an acetic acid-acidified, sodium 
chlorite aqueous solution at 70°C and then rinsed thor- 
oughly in DI water. Final isolation of a-cellulose was 
achieved by treatment in 17% NaOH solution following 
methodology identified in Sternberg (1989) with samples 
being rinsed in DI water and dried at 70°C. 

The a-cellulose was combusted to C0 2 in a Thermo- 
Finnigan TC/EA and then delivered to a Thermo Finnigan 
Delta Plus XL mass spectrometer (Thermo Fisher Scien- 
tific Inc., Waltham, MA) operating in a flow-through 
mode with a Conflo III interface. We determined isotopic 
composition with respect to the PDB standard (Coplen 
1996). For every four to five samples, a homogenous acet- 
anilide working standard of known isotopic composition 
was run for mass spectrometer calibration to insure that 
machine precision of ca. 0.02-0.1 1% 0 (SD) was present 
for each batch of 30^0 samples analyzed at the Univer- 
sity of Arizona, Laboratory of Tree-Ring Research. Addi- 
tionally, we used a holocellulose laboratory standard 
(MAWS) for approximately every 15 samples, which 
established precision of 0.12% o . 

We determined intrinsic water-use efficiency (iWUE) 
following the equation of Ehleringer and Cerling (1995): 

iWUE = A/g = ( Ci /c a )0.625, (1) 

where A = C0 2 assimilation of ponderosa pine leaves, 
g = leaf stomatal conductance, 0.625 = constant, and cj 
c a = intercellular C0 2 and atmospheric C0 2 . We interpo- 
lated annual c a with decadal observations beginning in 
1800 (Etheridge et al. 1998) and with measurements col- 
lected at Mauna Loa from 1959 to 2009 (Keeling et al. 
2008). Equations 2 (carbon isotopic discrimination for 
plants A) and 3 (the q/c a ratio; Farquhar et al. 1982) 
allowed us to obtain (q): 

A = (<5 13 C a - <5 13 C P )/(1 + r5 13 C p /1000) (2) 

a/c a = {A-a)/{b-a) (3) 

We derived the stable carbon isotopic composition of 
the atmosphere (<5 13 C a ) from 1796 to 2002 using data 
from Francey et al. (1999) and Allison et al. (2003). 
Because yearly data were temporally incomplete until 
1991, we used a sixth-order polynomial to create annual 
resolution data (e.g., Hemming et al. 1998). Holding all 
pre- 1850 data stable at -6.4% 0 , we averaged <5 13 C a values 
by pentad. We also averaged isotopic composition of the 
trees that were pooled by site (8 13 C p ) by pentad. Respec- 
tively, a and b represented constants for the discrimina- 
tion during diffusion of C0 2 into the air (4.4% D ) and 
discrimination during carboxylation (assumed at 27% 0 ). 
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Statistical analysis 

We used Spearman correlation (PASW Statistics 18.0) to 
identify the principal climatic drivers of radial growth for 
each site paired with data from climatic divisions (NCDC 
1994) Montana- 1 (MT1, Kitchen Gulch) and Montana-4 
(MT4, Cabin Gulch) during 1895-2007/8. We examined 
monthly and seasonal precipitation, temperature, and 
drought severity influences, including lagged influences 
up to 1 year. 

We examined potential possibilities for increased MPB 
susceptibility by comparing BAI values between visually 
healthy and MPB-infected populations at each site during 
stress and non-stress years using the non-parametric 
Mann-Whitney U- test. During 1895-2007/8, we catego- 
rized all years with Palmer Drought Severity Index (PDSI) 
values <— 2.0 (moderate drought or greater) as stressful, 
or PDSI values > —0.5 (normal to wet) as non-stressful. 
We identified potential shifts in iWUE at Cabin Gulch 
following the methodology of Rodionov (2004) and Ro- 
dionov and Overland (2005). We used a significance level 
of P < 0.05 and cut-off length of 30 years to detect possi- 
ble multi-decadal-length changes in iWUE since 1860. 

Results 

Climate/growth relationships 

Climate-growth relationships indicated that the principal 
climate driver of growth was significantly correlated 
(P < 0.05) with July PDSI for each chronology. At Cabin 
Gulch, there was a large difference between chronologies 
(r s = 0.42 CGR, n = 20 cores and r s = 0.62 CGI, n = 18), 
but less difference existed between chronologies at Kitchen 
Gulch (r s = 0.55 KGR, n = 39; r s = 0.53 KGI, n = 30). 
During 1895-2008, July PDSI for MT4 had a downward 
linear trend (r s = -0.34, P < 0.01) while July PDSI for 
MT1 exhibited no significant change (r s = —0.11 
P = 0.24), indicating increasingly drier summer condi- 
tions at Cabin Gulch and stable conditions at Kitchen 
Gulch. During the last decade (1999-2008), however, 
there was severe summer dryness at Cabin Gulch, with all 
years having negative PDSI values and a mean PDSI value 
of —2.61. Kitchen Gulch also experienced dry July PDSI 
conditions over the same period (x = —1.85), but 2 years 
had PSDI values > 0. 

Basal area increment 

Basal area increment for MPB-infected and visually healthy 
trees at each site closely matched during the majority of 
the 150-year period (Figs. 2, 3), but diverged during the 
last two-three decades. Examining differences using 
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Figure 2. BAI values (cm 2 yr 1 ) for the two chronologies developed 
at Cabin Gulch during 1860-2007. 

30-year intervals, mean BAI for visually healthy trees at 
Cabin Gulch was significantly less during 1890-1919 and 
significantly greater during 1980-2007 when compared 
with the MPB-infected sample (Table 2). At Kitchen 
Gulch, mean BAI was significantly greater during 1980- 
2008 when compared with the MPB-infected sample 
(Table 2) and for the average of two visually healthy 
chronologies compared with the average of the MPB- 
infected samples. The inflection point for radial growth 
divergence occurred ca. 1980 for Cabin Gulch, with the 
visually healthy trees maintaining higher BAI values for 
every year and the largest differences beginning in the 
mid-1990s (Fig. 2). A similar but less pronounced pattern 
occurred at Kitchen Gulch, with consistently higher BAI 
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Figure 3. BAI values (cm 2 yr -1 ) for the two chronologies developed 
at Kitchen Gulch during 1860-2008. 



1144 



© 201 3 The Authors. Ecology and Evolution published by John Wiley & Sons ttd. 



P. A. Knapp et al. 



Mountain Pine Beetle Selectivity in Ponderosa Pine Forests 



Table 2. BAI comparison between chronologies for MPB-infected and visually healthy trees at both sites and the two sites combined for 30-year 
periods (AD 1860-2007/8). Shaded boxes indicate significance at p < 0.05. 



Period (AD) Chronology Segment length (Years) Mean BAI Std. deviation (BAI) Independent samples M-W U-test (one-tailed) 



1860-1889 KGR 


30 




36.60 


8.77 




KGI 


30 




33.76 


8.80 


0.092 


1890-1919 KGR 


30 




36.84 


8.44 




KGI 


30 




34.81 


8.59 


0.145 


1920-1949 KGR 


30 




1 9.83 


8.41 




KGI 


30 




19.06 


7.92 


0.351 


1950-1979 KGR 


30 




26.33 


7.90 




KGI 


30 




24.42 


7.63 


0.1 59 


1980-2008 KGR 


29 




19.37 


5.09 




KGI 


29 




16.62 


4.12 


0.025 


1860-1889 CGR 


30 




11.08 


2.89 




CGI 


30 




12.51 


3.81 


0.07 


1890-1919 CGR 


30 




9.43 


3.18 




CGI 


30 




11.59 


4.47 


0.02 


1920-1949 CGR 


30 




8.86 


2.03 




CGI 


30 




9.36 


2.50 


0.21 


1950-1979 CGR 


30 




10.24 


2.29 




CGI 


30 




9.37 


2.59 


0.12 


1980-2007 CGR 


28 




10.45 


2.87 




CGI 


28 




7.82 


2.67 


0.001 


Sites combined 


1 860-1 889 


Visually healthy 




30 


23.82 


5.55 




Infected 




30 


23.14 


5.89 0.282 


1890-1919 


Visually healthy 




30 


23.13 


5.42 




Infected 




30 


23.20 


6.14 0.477 


1 920-1 949 


Visually healthy 




30 


14.34 


4.71 




Infected 




30 


14.21 


4.59 0.474 


1950-1979 


Visually healthy 




30 


18.28 


4.69 




Infected 




30 


16.90 


4.63 0.131 


1980-2007 


Visually healthy 




28 


15.06 


3.46 




Infected 




28 


12.29 


2.97 0.002 


values for the visually healthy trees beginning 


5 in 


the early 


infected samples 


during 1860-2009 (101.5 /xmol/mol vs. 


1990s and the largest differences occurring durin; 


; the last 


100.9 /imol/mol) and 1980-2009 (118.3 vs. 117.7). Both 


decade (Fig. 3). The 


combined chronologies showed no 


samples experienced upward trends in iWUE, with signifi- 


differences between infected and healthy trees until 1980- 


cant (P < 0.05, 


30-year cut-off) shifts toward higher 


2007. 








iWUE beginning 


in 1955-59 for the visually healthy trees 


Mean BAI values 


at Cabin Gulch were 


18% higher 


and 1960-1964 


for the infected trees (Fig. 4). Mean 


(n = 27 years; P < 0.01, one-tailed) for 


the 


visually 


iWUE between the early and shift points was 96.3 vs. 


healthy trees compared with MPB-infected trees during 


113.75 (18% increase) and 96.2 vs. 110 (17% increase) for 


stressful years (PDSI< 


—2.0), but the two populations were 


the visually healthy and the MPB-infected populations. 


not different (n = 66; 


P = 0.781, two-tailed) during non- 






drought years (PDSI > -0.5). At Kitchen Gulch, BAI val- 
ues of healthy trees were not significantly greater during 


Discussion 




either drought (n = 
non-drought years (n 


26 years, P = 0.19, one-tailed) or 
= 69; P = 0.17, two-tailed). 


Radial growth divergence 










Basal area increment values within populations and 


Intrinsic water-use efficiency 






between sites were similar until the late 20 th Century, at 
which point the MPB-infected populations had consis- 


Intrinsic water-use efficiency at Cabin Gulch 


was 


not sig- 


tently lower BAI values. Thus, growth rates 20-30 years 


nificantly different (P > 0.05) for visually healthy and 


prior to the current outbreak diverged between our 
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Figure 4. Intrinsic water-use efficiency (umol/mol) for visually healthy 
(CGR) and MPB-infected trees (CGI) at Cabin Gulch (pentads from AD 
1860-2009). Vertical arrows indicate regime shift (p < 0.05) inflection 
points. 

selected populations with the slower growing trees being 
more vulnerable to infestation. These findings are consis- 
tent with the vigor hypothesis (e.g., Larsson et al. 1983), 
and the findings may be revealing in that these were 
open-grown trees sampled from the same slope, soils, and 
aspect as visually healthy trees. Sampled trees from the 
chronologies were often (approximately 50% of total) 
collected within 15 m of each other, further reducing the 
likelihood that growth differences were an artifact of 
microsite variability. 

Although changes to higher iWUE occurred slightly 
earlier for the visually healthy chronology, iWUE trended 
upward for both chronologies at Cabin Gulch with only 
minor and non-significant differences between popula- 
tions. Thus, changes in iWUE do not explain the growth 
divergences. Growth/climate responses also were similar, 
with both populations at both sites dependent on sum- 
mer soil-moisture conditions as inferred from July PDSI 
values. At Cabin Gulch, the infected population exhibited 
greater growth sensitivity to climate, which is consistent 
with the findings of McDowell et al. 2010, who found 
greater drought-induced mortality in ponderosa pine with 
increasing climatic sensitivity. Conversely, at Kitchen 
Gulch, the healthy population exhibited slightly greater 
sensitivity to climatic variations than the infected popula- 
tion, suggesting that the relationship between vigor and 
climatic sensitivity may be site specific. 

Intrinsic water-use efficiency 

Our findings also suggest that significant increases in 
iWUE do not confer immunity to MPB infestations for 



all trees during moisture-stress periods, but may have the 
potential to lessen the severity of the attacks because of 
the ameliorative effects of improved plant-water relation- 
ships. In the Northern Rockies, USA, radial growth rates 
and iWUE of ponderosa pine are significantly and 
positively associated with age, suggesting that the water 
stress-ameliorating effects of increased atmospheric C0 2 
preferentially favor older trees (Knapp and Soule' 2011). In 
this study, median tree age at CGR (260 years) and KGR 
(180 years) were greater than for CGI (207 years) and 
KGI (156 years). Our findings are consistent in that older 
trees may have been less affected by drought conditions 
than younger trees during the last 30 years, and in turn 
have been less susceptible to MPB infestations due to the 
reduction in moisture stress, particularly at Cabin Gulch 
where age differences (53 years) were more pronounced. 

Implications 

The divergence in growth shown between populations at 
these sites may reveal "previously invisible genetic quirks" 
(Neese and Williams 1994:144) that may favor/disfavor 
different genotypes embedded within a local population 
when stress occurs. For both sites sampled, MPB-infected 
and visually healthy populations followed two different 
growth trajectories in the late 20th Century/early 21st 
Century, suggesting that some type of environmental shift 
favored one subgroup of ponderosa pine over the other 
prior to the onset of MPB attacks. In turn, the less-vigor- 
ous trees were unable to repel MPB attacks, making them 
more vulnerable to the current epidemic. For the Cabin 
Gulch and Kitchen Gulch chronologies, the visually 
healthy trees were less affected by drought conditions as 
they maintained higher BAI values during the past 20- 
30 years, although this relationship was better expressed 
at Cabin Gulch where the overall climatic conditions are 
drier than at Kitchen Gulch and age differences between 
groups were less. Conversely, the lack of pronounced 
superior growth of the healthy trees throughout the 
record suggests other possibilities for divergence. For 
example, prior MPB infestations that weakened, but did 
not kill, the trees 20-30 years ago may have made the 
MPB trees more susceptible to the current episode. 

Overall these results were similar between the Cabin 
Gulch and Kitchen Gulch sites. Radial growth divergence 
occurred several decades prior to infection, suggesting 
that the likelihood of MPB infestation of individual trees 
older than 150 years is not arbitrary and potentially set in 
motion decades prior to actual infestations. In short, sus- 
ceptibility of ponderosa pine to MPB infestations may be 
strongly associated with the long-term vigor of individual 
trees that emerges during episodic outbreaks. The conse- 
quences of these outbreaks are that they selectively 
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remove the least vigorous trees from the local forest pop- 
ulation, thus promoting an overall increase in tree vigor 
and drought resistance. 
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